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Bes chr eibung 

Verfahren und Anordnung sowie Computerprograinni mit Programm- 
code-Mitteln und Computerprograinm-Produkt zur Ermittlung ei- 
5 ner ausgewahlten Position einer mobilen Koinniunikationsein- 
richtung in einem Koznmunikationsnetz 

Die Erfindung betrifft eine Ermittlung einer Position einer 
mobilen Konimunikationseinrichtung in einem Kommunikationsnetz 
10 (Lokalisierung) • 

Mit zunehmender Verbreitung mobiler Kommunikation nimmt auch 
eine Nachfrage nach zusatzlichen Dienstleistungen, sogenann- 
ten "^services", bei Mobil funksystemen zu. 

15 

Unter "^Location Based Services'' werden dabei zusatzliche 
Dienstleistungen von Mobilfunkbetreibern verstanden, welche 
Nutzern von den Mobilfvinksystemen ortsgebunden^ d.h. abhSngig 
von einer Position oder einem Standort des jeweiligen Nut- 
20 zers, angeboten bzw. bereitgestellt werden, beispielsweise 

orts- Oder entf ernungsabhangige Nutzungstarif e Oder Orientie- 
rungshilfen fUr Rettungseinsatze oder Suchdienste. 

Grundlegend fUr einen "^Location Based Service" ist demzufolge 
25 eine Lokalisierung oder Positionsermittlung des jeweiligen 
Nutzers bzw. dessen mobilen Kommunikationsgerats- 

Ftlr eine solche Lokalisierung von mobilen Kommunikationsgera- 
ten in Kommunikationsnetze sind verschiedene Techniken be- 

30 kannt, beispielsweise eine Positionsermittlung auf Basis ei- 
ner Lauf zeitbestimmung bzw. Lauf zeitmessung von Kommunikati- 
onssignalen von einem mobilen Kommunikationsgerat zu einer 
Basisstation eines Kommunikationsnetzes ([1], [2]) oder eine 
Lokalisierung mittels satellitengesttitzter Systeme, wie ein 

35 GPS- 
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Das aus [2] bekannte, Laufzeit basierte Positionsenaittliings- 
verfahren wird fiir ein Mobiltelef on, im allgemeinen einer Mo- 
bilstation, in einem GSM-Koinmunikationsnetz (= Global System 
for Mobile Communications) ([7], [8], [9]) gemafl einer TDMA- 
5 Mobil funk-Technik (Time Division Multiple Access-Mobilfunk) 
durchgef tihrt . 

Einer einzelnen Mobilstation, die sich bei einer ortsfesten 
Basisstation (gesprachsftihrenden Basisstation) eingebucht 
10 hat/ wird bei dieser eine derzeit freie Zeitlage in einem 
TDMA-Rahmen zugeteilt. 

Zu dieser Zeitlage gehen die ftlr die betreffende Mobilstation 
bestimmten Kommunikationssignale in Signalpaketen, sogenann- 
15 ten bursts, mit einer Lange von 15/26ms von der Basisstation 
ab, beziehungsweise mUssen die von der Mobilstation gesende- 
ten Kommunikationssignale bzw. bursts an der Basisstation 
eintreffen. 

20 Die von der Basisstation emittierten Kommunikationssignale 

finden ihren Weg zur Mobilstation aufgrund von Streuungen tim- 
ber unterschiedliche Wege (Mehrf achausbreitungen) , wobei sie 
frequenzabhangig gedampft werden. 

25 Eine Empf angsf eldstarke der von der Mobilstation empfangenen 
Kommunikationssignale ist somit nicht nur von einer Entfer- 
nung der Mobilstation von der Basisstation, sondern auch von 
der Frequenz xand topographischen. Gegebenheiten zwischen Mo- 
bilstation und Basisstation abhangig. Deshalb werden die ein- 

30 zelnen Datenpakete auf verschiedenen Tragerf recpienzen ver- 

schickt, wodurch selektive Storungen einer Frequenz auf meh- 
rere Teilnehmer verteilt werden kSnnen. 

Daftir ist allerdings eine prazise Synchronisierung zwischen 
35 Mobilstation und Basisstation notwendig. Diese Synchronisie- 
rung wird durch eine Mobilitat eines Nutzers zusatzlich er- 
schwert, well die Mobilstation sich nun in sich andernden 
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Entfernungen zur Basisstation befindet und ihre Kommunikati- 
onssignale unterschiedliche Laufzeiten aufweisen. 

Urn die unterschiedlichen Laufzeiten auszugleichen und der Ba- 
5 sisstation rahmensynchrone Daten liefern zu kSnnen/ misst die 
Mobilstation die Signallauf zeit zur Basisstation und korri- 
giert dadurch einen Sendebeginn ihrer bursts. 

Die Signallauf zeit wird in einem sogenannten timing advance" 
10 (TA) codiert und weist eine Abhangigkeit zur Entfernung zwi- 
schen Mobilstation und gesprachsftlhrenden Basisstation auf . 

FUr das TA stehen 64 Stufen zur Verftigung^ die mit Werten 0 
bis 63 (bit-) codiert werden und die Laufzeit reprasentieren. 

15 

Da Positionen von Basisstationen bekannt sind, kann aus einem 
TA bzw. aus der Signallauf zeit auf die Position der Mobilsta- 
tion geschlossen werden. 

20 Eine Messgenauigkeit bei der Lauf zeitbestimmung betrSgt eine 
Bitdauer, also in GSM 48/13 \is, was einer einfachen Wegiange 
von etwa 554m entspricht. 

Aus [3] ist eine Positionsermittlung einer mobilen Kommunika- 
25 tionseinrichtung in einem UMTS-Netz (= Universal Mobile Tele- 
communication System-Netz) bekannt . 

Bei dem entsprechenden UMTS- Mobil funks tandard, auf welchem 
das UMTS-Netz basiert/ ist eine Positionsermittlung eines Mo- 
30 bilfunkgerates bereits explizit im Standard enthalten bzw. 
wird durch diesen gefordert (TS 25.305 V3.1.0: stage 2 
^'Functional Specification of Location Services in UTRAN" (re- 
lease 99), 3GPP TSG-RAN-WG2, 2000). 

35 Aus [4], [5] und [6] sind weitere Verfahren zur Lokalisierung 
einer mobilen Kommunikationseinrichtung in einem Kommunikati- 
onsnetz bekannt. 
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Aus [10] ist ein nichtlineares Filter, ein Prior Density 
Splitting Mixture Estimator (PDSME) / bekannt. 

5 Dieses nichtlineare Filter, dieser PDSME, basiert auf neue 
Gaussche Mischfilter-Algorithmen fUr nichtlineare Systemzu- 
standsgleichungen und/oder nichtlineare Messgleichungen. 

Oblicherweise treten Linearisierungsf ehler auf, wenn die 
10 nichtlinearen Systemzustandsgleichungen und/oder die nichtli- 
nearen Messgleichungen durch Linearisierungen ersetzt werden 
bei Anwendung von linearen Filtertechniken. 

Der PDSME bewirkt eine Reduktion dieser Linearisierungsf ehler 
15 durch Splitten der a priori Dichte in mehrere Gaussche Misch- 
komponenten mit jeweils kleineren Kovarianzen. 

Der PDSME kann angewandt werden sowohl fUr einen sogenannten 
Prediction Step" bei den nichtlinearen Systemzustandsglei-- 
20 chungen als auch ftir einen sogenannten Filter Step'' bei den 
nichtlinearen Messgleichungen. 

Weiter ist aus [10] bekannt eine Messung ftlr den Linearisie- 
rungsfehler, die vergleichbar einer Kullback-Leibler-Distanz 
25 ist. Dartlber kann eine gewUnschte Abschatzungsgtite bei dem 

PDSME eingestellt werden. Durch BeschrSnkung einer maxiiaalen 
Anzahl von Mischkomponenten bei dem PDMSE kann eine obere 
Grenze ftir einen Rechenaufwand angegeben werden. 

30 Die genannten Lokalisierungsverf ahren weisen unter andereia 

die Nachteile auf, dass die durch sie ermittelten Positionen 
der mobilen Kommunikationseinrichtungen ungenau und demzufol- 
ge mit groiien Unsicherheiten behaftet sind. Genauere Verf ah- 
ren erfordern allerdings aufwendige Zusatzeinrichtungen und 

35 kostspielige Modif ikationen an den bzw. der Kommunikations- 
netze und Kommunikationseinrichtungen. 



wo 2004/095868 



5 



PCT/EP2004/003320 



Somit liegt der Erfindung die Aufgabe zugrunde^ eine genaue 
und mit moglichst geringen Unsicherheiten behaftete Lokali- 
sieriing einer mobilen Koitununikationseinrichtung in einem Kom- 
munikationsnetz zu ermoglichen/ welche moglichst einfach und 
5 kostengtinstig realisiert werden kann. 

Diese Aufgaben werden durch das Verfahren und die Anordnung 
sowie durch das Compute rprogramm lait Prograitaacode-Mitteln und 
das Computerprogramm-Produkt zur Ermittlung einer ausgewahl- 
10 ten Position einer mobilen Kommunikationseinrichtung in einem 
Kommunikationsnetz mit den Merkmalen gemaB dem jeweiligen un- 
abhcingigen Patentanspruch gelSst. 

Bei dem Verfahren zur Eirmittlung einer ausgewahlten Position 
15 einer mobilen Kommunikationseinrichtung in einem Kommunikati- 
onsnetz mit mindestens einer ersten Basisstation, eingerich- 
tet zu einer ersten Kommunikation mit der mobilen Kommunika- 
tionseinrichtung mittels erster Koromunikationssignale, wird 
ein nichtlineares Kommunikationsmodell mit einem determinis- 
20 tischen Anteil und einem stochastischen Anteil verwendet. 

Dieses nichtlineare Kommunikationsmodell umfasst einen deter- 
ministischen und einen stochastischen Anteil 

- welcher deterministische Anteil eine AbhSngigkeit zwischen 
25 Kommunikationssignale einer Basisstation \ind Positionen 

der mobilen Kommunikationseinrichtung beschreibt und 

- welcher stochastische Anteil eine Unsicherheit des deter- 
ministischen Anteils beschreibt. 

30 Unter Verwendung des zu der ausgewahlten Position zugehSrigen 
ersten Kommunikationssignals der ersten Basisstation und dem 
nichtlinearen Kommunikationsmodell wird die ausgewahlte Posi- 
tion ermittelt, wobei bei der Ermittlung eine nichtlineare 
Bayesianische Filtertechnik angewandt wird- 

35 

Durch Anwendung der nichtlinearen Bayesianische Filtertechnik 
auf das Kommunikationsmodell bzw. auf das Kommunikationssig- 
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nal wird ein mogliches Auf enthaltsgebiet der ausgewahlten Po- 
sition bestimmt, aus welchem schlieBlich die ausgewahlte Po- 
sition bestimmt wird. 

5 Unter der nichtlinearen Bayesianischen Filtertechnik bei der 
Erfindtmg ist im allgemeinen die Technik wie folgt zu verste- 
hen: 

- die Bayesianische Filtertechnik entspricht einer zeitli- 
chen Fortschreibung bzw. Veranderung (auch Berechnung der 
10 zeitlichen Veranderung) bekannter Wahrscheinlichkeitsver- 

teilungen der SystemzustSnde und Systemausgange durch das 
Vorliegen neuer Messwerte. 

Dem erf inderischen Lokalisierungsverf ahren liegt die Idee 
15 zugrunde, aus in Kommunikationsnetzen, wie zum Beispiel WIAN, 
GSM, DECT-Netzen, zur VerfUgung stehenden Koinmunikationssig- 
nalen zwischen mindestens einer Basisstation (Sender) und ei- 
ner Mobilstation (Empfanger) entf ernungsrelevante Parameter 
und daraus geographische Inf ormationen, in diesem Fall ein 
20 mSgliches Aufenthalts- oder Entfernungsgebiet bzw. Unsicher- 
heitsgebiet der Mobilstation, zu gewinnen. 

Anzumerken ist, dass Sender und Empfanger der Kommunikations- 
signale technisch auch vertauscht realisiert sein kOnnen. Die 
25 erfinderische Vorgehensweise bleibt davon unbertihrt. In die- 
sem Fall ware die Basisstation der Empfanger und die Mobil- 
station der Sender der in Rede stehenden Kommunikat ions signa- 
le. 

30 Weiter liegt der dem erf inderischen Lokalisierungsverf ahren 
die Erkenntnis zugrunde, dass die von der Mobilstation emp- 
fangenen Kommunikationssignale abhangig sind von der geogra- 
phischen Position der Mobilstation im Kommunikationsnetz bzw. 
zu der Basisstation. 

35 

Diese AbhSngigkeit zwischen den empfangenen Kommunikations- 
signalen bzw. zwischen Signalstarken der (von der Basisstati- 
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on) empfangenen Kommunikationssignale land der geographischen 
Position der Mobilstation kann durch sogenaxinte Messgleichun- 
gen, d.h, durch das Kommunikationsmodell^ beschrieben werden. 

5 Diese Messgleichungen bzw. das Kommunikationsmodell umfassen 
einen deterministischen Anteil, welcher deterministische An- 
tell eine Abhangigkeit zwischen dem empfangenen Kommunikati- 
onssignal nnd der Position der Mobilstation beschreibt# und 
einen stochastischen Anteil, welcher stochastische Anteil ei- 
10 ne Unsicherheit des deterministischen Anteils beschreibt. 

Diese Unsicherheit kann beispielsweise eine solche des Kommu- 
nikationssignals und/oder eine solche der erwShnten Abhangig- 
keit (Modellunsicherheit) sein, 

15 

Die erf indungsgemafi auf diese Messgleichungen bzw. das Kommu- 
nikationsmodell angewandte nichtlineare Bayesianische Filter- 
technik zur Bestimmung des m5glichen Auf enthaltsgebiets der 
Mobilstation bzw. der ausgewahlten Position verbessert die 
20 Qualitat der Lokalisierung verglichen mit Ublichen Filter- 

techniken wie einem Extended Kalmann Filter (EKF) erheblich. 

Die Anwendung der nichtlinearen Bayesianischen Filtertechnik 
kann in einem sogenannten Filter Step" einmalig oder itera- 
25 tiv mehrmalig erf ol gen. 

In diesem mOglichen Auf enthaltsgebiet bzw. Unsicherheitsge- 
biet/ beschrieben unter Verwendung stochastischer GrdBen, wie 
einem Mittelwert und/oder einer Varianz, wird schliefilich die 
30 Mobilstation bzw. deren Position vermutet. 

Diese Vermutung kann mathematisch formuliert werden durch ei- 
nen Kennwert des moglichen Aufenthalts- bzw. Unsicherheitsge- 
biets, wie einem Schwerpunkt oder einem Erwartungswert . Die- 
35 ser kann dann als Schatzung ftir die ausgewahlte Position der 
Mobilstation verwendet wird. 
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Ein besonderer Vorteil der Erfindung liegt darin, dass die 
Lokalisierung auf Basis von Koininunikationssignalen \ind be- 
kannten Positioner! von Basisstationen durchgeftihrt wird, wel- 
che in einem nonaalen Betrieb bei einem Mobilfunksystem an- 
5 fallen und dort zur Verftigung stehen. Dadurch* kann auf auf- 
wendige Anderungen und Erweiterungen sowie Zusatzmessungen 
bestehender Mobilfunksysteme bzw. bei bestehenden Mobilfunk- 
systemen verzichtet werden. 

10 Das Computerprogramm mit Programmcode-Mitteln ist eingerich- 
tet, urn alle Schritte gemaB dem erf indungsgemafien Verfahren 
zur Ermittlung einer Position, d,h, dem erf indungsgemSBen Lo- 
kalisierungsverfahren, durchzuftihren, wenn das Programm auf 
einem Computer ausgefUhrt wird. 

15 

Das Computerprogramm-Produkt mit auf einem maschinenlesbaren 
Trager gespeicherten Programmcode-Mitteln ist eingerichtet, 
urn alle Schritte gemai3 dem erf indungsgemSBen Lokalisierungs- 
verfahren durchzuftihren, wenn das Programm auf einem Computer 
20 ausgeftihrt wird. 

Die Anordnung sowie das Computerprogramm mit Programmcode- 
Mitteln, eingerichtet urn alle Schritte gemSfl dem erfinderi- 
schen Lokalisierungsverf ahren durchzuftihren, wenn das Pro- 

25 gramm auf einem Computer ausgeftihrt wird, sowie das Computer- 
pro gramm-Produkt mit auf einem maschinenlesbaren Trager ge- 
speicherten Programmcode-Mitteln, eingerichtet um alle 
Schritte gemaB dem erf inderischen Lokalisierungsverf ahren 
durchzuftlhren, wenn das Programm auf einem Computer ausge- 

30 ftihrt wird, sind insbesondere geeignet zur Durchftihrung des 
erf indungsgemaBen Lokalisierungsverf ahrens Oder einer seiner 
nachfolgend eriauterten Weiterbildungen. 

Bevorzugte Weiterbildungen der Erfindung ergeben sich aus den 
35 abhangigen Ansprtichen. 
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Die im weiteren beschriebenen Weiterbildungen beziehen sich 
sowohl auf die Verfahren als auch auf die Anordnung. 

Die Erfindung und die im weiteren beschriebenen Weiterbildun- 
5 gen konnen sowohl in Software als auch in Hardware/ bei- 
spielsweise iinter Verwendung einer speziellen elektrischen 
Schaltung^ realtsiert werden. 

Ferner ist eine Realisierung der Erfindiing oder einer im wei- 
10 teren beschriebenen Weiterbildung mOglich durch ein computer- 
lesbares Speichermediiim/ auf welchem das Computerprogramm mit 
Programmcode-Mitteln gespeichert ist/ welches die Erfindung 
Oder Weiterbildung ausftlhrt. 

15 Auch kann die Erfindung oder jede im weiteren beschriebene 

Weiterbildung durch ein Computerprogrammerzeugnis realisiert 
sein/ welches ein Speichermediiam aufweist/ auf welchem das 
Computerprogramm mit Programmcode-Mitteln gespeichert ist/ 
welches die Erfindung oder Weiterbildung ausfUhrt, 

20 

Bei einer Kommunikation in einem Kommunikationsnetz zwischen 
einer mobilen Kommunikationseinrichtung (Mobilstation) , bei- 
spielsweise einem Mobiltelefon, und einer Basisstation/ bei- 
spielsweise einer Rundantenne bzw. einem Rundstrahler oder 
25 einer bzw. mehrerer sektoraler Antennen/ werden Daten/ die 
Kommuni kat ions signal e/ in Signalpaketen/ sogenannten bursts, 
tibertragen. 

Basierend auf den bzw. unter Verwendung der tlbertragenen Kom- 
30 munikationssignalen bzw. Signalpaketen lassen sich verschie- 
dene entf ernungsrelevante Parameter ermitteln, welche wieder- 
um als Grundlage fUr die Ermittlung von den moglichen Aufent-- 
halts- bzw. Entfernungsbebieten herangezogen werden konnen. 

35 Ein solcher entf ernungsrelevanter/ d.h. entf ernungsabhangi- 
ger/ Parameter ist beispielsweise eine Feldstarke eines Sig- 
nalpakets . 
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Die Feldstarke weist eine natiirliche Abhangigkeit zur Entfer- 
nung zwischen der Mobilstation und der (gesprSchsfiihrenden) 
Basisstation auf und liefert demzufolge eine Information iiber 
5 das m5gliche Auf enthaltsgebiet bzw. Entf ernungsgebiet (Unsi- 
cherheitsgebiet) der Mobilstation. 

Diese AbhSngigkeit zwischen Feldstarke und Entfernung kann 
durch physikalische Modelle, welche ein Ausbreitungsverhalten 
10 von Signalen beschreiben, beschrieben werden, 

FUr die Ermittlung der ausgewahlten Position kann ein Kenn- 
wert des moglichen Aufenthalts- bzw. Unsicherheitsgebiets be- 
stimmt werden, wie ein Schwerpunkt Oder ein Erwartungswert, 
15 welcher dann als Schatzung fUr die ausgewahlte Position der 
Mobilstation verwendet wird. 

Die Qualitat der Lokalisierung kann welter verbessert werden 
unter Verwendung eines Benutzermodells, welches eine Bewegung 
20 der mobilen Kommunikationseinrichtung beschreibt. So kann ein 
solches Benutzermodell eine maximale Schrittweite bzw. Bewe- 
gungsdistanz ftlr einen vorgegebenen Zeitschritt begrenzen. 

Die Erfindung eigne t sich insbesondere zu einem Einsatz im 
25 Umfeld eines digit alen, zellularen Mobilfunksystems/ wie ei- 
nes GSM-Netzes/ und dort beispielsweise zur Lokalisierung ei- 
nes GSM— Telefons (Mobiltelefon) . 

Dabei werden bei dem Einsatz der Erfindung nur die dem Mobil- 
30 telefon zur Verfiigung stehenden Daten verwenden, wobei weder 
am GSM-Netz noch an Mobilstationen in dem GSM-Netz kostspie- 
lige Anderungen vorzunehmen sind. 

Beispielsweise sind bekannt von einem GSM-Netz die Positionen 
35 der einzelnen Basisstationen und deren Antennen sowie deren 
Charakteristika, welche Auskunft iiber das Versorgungsgebiet 
der jeweiligen Antenne gibt. 
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Das Mobiltelefon seinerseits steht far einen korrekten Ver- 
bindvingsaufbau immer in Kontakt mit den 

empf angbaren Antennen, urn die am besten fUr ein Gesprach 
5 geeignete Antenne vom Netz zugeteilt zu bekoiniaen. Dazu misst 
es unter anderem die Empf angsf eldstarke der empfangbaren Rn- 
tennen sowie bestimmt Signallauf zeiten, die dann ebenfalls 
bekannt sind. 

10 Auf Basis dieser vorliegenden Informationen erfolgt dann die 
Lokalisierung des Mobiltelefons gemaB der erf inderischen Vor- 
gehensweise. 

Auch eignet sich die Erfindung zu einem Einsatz im Umfeld 
15 weiterer digitaler, zellularen Mobilfunksysteme/ wie eines 

WLAN Oder eines DECT-Netzes [11] f und dort beispielsweise zur 
Lokalisierung eines DECTHMobiltelefons . 

In Figuren ist ein Ausftihrungsbei spiel der Erfindung darge- 
20 stellt/ welches im weiteren naher eriautert wird. Das Ausftlh- 
rungsbeispiel ist dabei unterteilt in einen grundlegenden 
Teil, welcher Grundlagen der erf inderischen Vorgehensweise 
darstelltr und einen angewandten Teilr welcher anhand eines 
konkreten Zahlenbeispiels die erf inderische Vorgehensweise 
25 sowie Ergebnisse verdeutlicht bzw. angibt. 

Es zeigen 

Figur 1 eine Skizze, in welcher eine angenommene lineare NS- 
30 herung ftir das nichtlineare Kommunikationsmodell dar- 

gestellt ist; 

Figur 2 eine Skizze, in welcher die Unsicherheit der Kommuni- 
kationssignale durch rauschen dargestellt ist; 

35 

Figur 3 eine Skizze, in welcher ein Oberblick tiber den PDSME- 
Algorithmus mit einem linearisierten Fiterstep (obe- 
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rer Ausschnitt) und einem Prognoseschritt (Predikti- 
onsstep) (unterer Ausschnitt) dargestellt ist; 

Figuren 4a und b Skizzen, in welchen ein Evaluierungs- 
5 Szenario fur die Lokalisierung gemafi dem AusfUhrungs- 

beispiel dargestellt ist; 

Figuren 5a vind b Skizzen, in welchen eine Lokalisierung unter 
Verwendung eines Extended Kalman Filters dargestellt 
10 ist (erster Filterschritt Fig. 5a; 170-ter Filter- 

schritt Fig, 5b) ; 

Figuren 6a und b Skizzen, in welchen eine Lokalisierung unter 
Verwendung des PDSME gemafi dem Ausftihr\ingsbeispiel 
15 dargestellt ist (erster Filterschritt Fig. 6a; 170- 

ter Filterschritt Fig. 6b; 

Figuren 7a und b Skizzen, in welchen Ergebnisse der EKF- 

Filterung und der PDSME-Filteriing gemafi dem Ausftlh- 
20 rungsbeispiel in Bezug auf die korrekten Ergebnisse 

dargestellt sind; 

Figur 8 eine Skizze, in welcher eine Anzahl von Gausschen 
Mischkomponenten in einer Abhangigkeit von Filter- 
25 schritten dargestellt sind. 



Ausftihrungsbel spiel : Lokalisierung eines DECT-Mobiltelefons 
in einem DECT-Netz basierend auf einer nichtllnearen Filter- 
30 technik 

Grundlagen 

Nachfolgend wird die Anwendung einer nichtllnearen Bayesiani- 
35 schen Filtertechnik bei der Lokalisierung eines mobilen Funk- 
telefons beschrieben. Die Anwendung wird beschrieben am Bei- 
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spiel einer Lokalisierung eines DECT-Mobiltelefons in einem 
DECT-Netz mit mehreren Basisstationen. 

Signalstarken (Feldstarken) der von einem Mobil telef on emp- 
5 fangenen und dort gemessenen Kommunikationssignale (einer 
diese Signale aussendenden Basisstation) sind abhangig von 
der Position des Mobiltelefons beztiglich der die Konimunikati- 
onssignale aussenden Basisstation. Dieser Zusammenhang wird 
beschrieben durch nicht lineare Messgleichungen (nichtlinea- 
10 res Kommunikationsmodell) . 

Diese nicht linearen Messgleichungen bzw. dieses nicht linea- 
re Kommunikationsmodell umfasst einen deterministischen An- 
teil/ welcher die empfangenen Signalstarken als Funktion der 
15 Position beschreibt, und einen stochastischen Anteil, welcher 
Modellfehler und Messrauschen berticksichtigt . 

ZusStzlich werden Benutzermodelle mit berUcksichtigt , welche 
ein Wissen iiber eine (raimliche) Bewegung von Benutzern der 
20 Mobiltelefone mit einbringen. 

Eine neue nichtlineare Filtertechnik, ein sogenannter Prior 
Density Splitting Mixture Estimator" (PDSME) , einordnbar als 
Gausscher Mischf ilteralgorithmuS/ verbessert signifikant eine 
25 Lokalisierungsqualitat verglichen mit Standardf iltertechni- 
ken, wie ein Extended Kalman FLltex" (EKF) (vgl. Figuren 5 
bis 7) . 

Typische Anwendungen fUr Lokalisierungsauf gaben sind Positi- 
30 onseinxiittlungen mobiler Kommunikationseinrichtungen (Mobilge- 
raten) in WLANs, GSM-Netzen oder DECT-Netz en. 

Ftir die Positionsermittlung sollten dabei Inf ormationen ver- 
wendbar sein, welche bereits wahrend einer normalen Operati- 
35 onsphase, d.h. bei einer "^normalen'' Kommunikation zwischen 
Mobilgerat und Basisstation, zur Verftigung stehen. 
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Das aus [12] bekannte Lokalisierungsverf ahren fUr Mobiltele- 
fone in GSM-Netzen basiert auf ortsabhangige/ von den Mobil- 
telefonen empfang- und dort messbare Signalstarken der von 
Basisstationen ausgesendeten Kommunikationssignale . Wegen der 
5 Ortsabhangigkeit sind diese charakteristisch fUr die Position 
des Signalempf angers, in diesem Fall ftlr die Position des Mo- 
biltelefons bzw. eines Benutzers des Mobiltelefons • 

wahrend der normalen Operationsphase werden die empfangenen 
Signalstarken aller empf angbarer Kommiinikationssignale der 
jeweiligen Sender (Basisstationen) durch das Mobiltelefon ge- 
messen xind verglichen, Dadurch wird ein sogenanntes Handover 
des Mobiltelefons zwischen verschiedenen Sendern bzw. Basis- 
stationen gesteuert. 

Basierend auf diese im Normalbetrieb anfallenden Messungen 
konnen Lokalisierungsverf ahren entwickelt werden, um zusatz- 
lich die Position des Mobiltelefons in Bezug auf die Sender 
bzw. im GSM-Netz zu ermitteln. 

Aus [13] und [17] sind weitere Vorgehensweisen zur Lokalisie- 
rung mobiler Kommunikationseinriclitungen bekannt. Diese Vor- 
gehensweise bei der Lokalisierung basieren auf einer Anwen- 
dungen von nichtlinearen Zustandsschatzungen. 

Nachfolgend wird die Anwendung einer neuen nichtlinearen Fil- 
tertechnik ftir die Lokalisierung von DECT-Mobiltelefonen be- 
schrieben. 

30 Um die Lokalisierung von Mobiltelefonen zu ermOglichen, wird, 
wie aus [18] bekannt ist, ein stochastisches Naherungsmodell 
fUr die Ausbreitung von Radiowellen Oder empfangbaren Signal- 
starken von Koromunikationssignalen in einer Lo kalis ierungsijun- 
gebung identif iziert bzw. entwickelt. 



35 
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Ein entsprechendes Ausbreitungsmodell kann auch grundsatzlich 
durch physikalische Ausbreitungsmodelle von elektromagneti- 
schen Wellen erlangt werden. 

5 Eine Anpassung rein physikalischer Ausbreitungsmodelle an re- 
ale (Lokalisierungs^>Umgebungen ist jedoch sehr kompliziert/ 
da Reflexionen oder Interferenzen oder andere physikalische 
Phanomene die physikalisch angenommene Ausbreitung beeinflus- 
sen und verfalschen. 

10 

Da reale Ausbreitungsparameter/ welche die elektromagneti- 
schen (Ausbreitung- )Eigenschaf ten in der realen Lokalisie- 
rungsumgebung, in welcher das Mobiltelefon zu lokalisieren 
ist/ beschreiben/ nur partiell bekannt bzw. bestimmbar sind/ 
15 ist die Anwendung solcher physikalisch basierter Ausbrei- 
tungsmodelle fiir solche realen Anwendungen ausgeschlossen. 

Deshalb ist ein Ausbreitungsmodell, welches auf tatsachlichen 
Messungen der (logarithmierten) empfangenen Signalstarken der 
20 Kommunikationssignale in der tatsachlichen Lokalisierungsiam- 
gebung beruht^ bzw. sich daraus ableiten lasst, theoretischen 
physikalischen Ausbreitungsmodellen vorzuziehen. 

25 Die Entwicklung eines solchen aus tatsachlichen Messungen ab- 
geleitetes ""Messmodell" bzw. Koinmunikation- bzw. Ausbrei- 
tungsmodell erf ordert ein Einmessen oder eine Kalibrierung im 
Rahmen dessen Parameter des Messmodells festgelegt werden. 

30 Eiiuaessen/Kalibrieren (""model generation phase") 

FUr ein solches Einmessen/Kalibrieren wird der Empfanger (be- 
nutzt im Sinne eines Messgerates) , d.h. das Mobiltelefon, an 
mehrere Positionen in der Lokalisierungsumgebung bewegt, bei- 
35 spielsweise an Gitterpunkten erhaltlich durch Platzierung ei- 
nes Gitters auf die Lokalisierungsumgebung. 
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An jedem Gitterpiinkt werden die Feldstarken aller dort emp- 
fangbaren Basisstationen gemessen und in einer Karte der Lo- 
kalisierungsumgebung gespeichert, Auf diese Weise wird die 
Feldstarkenverteilung jeder Basisstation tiber die Lokalisie- 
5 rungsiHagebung bestimmt. Die Feldstarkenverteilungen konnen 
durch das Messmodell modelliert werden. 

Basierend auf diese Modellinf ormationen (und von an einer ak- 
tuellen Position gemessenen Feldstarken) kann dann die aktu- 
10 elle/ von einem Empf anger bzw, von dem Mobil telef on in der 

Lokalisierungsumgebung eingenommenen Position bestimmt werden 
(Lokalisierung in einer ^localisation phase"). 

Weiter wird nachfolgend beschrieben werden^ wie ein sto- 
15 chastisches Messmodell in der model generation phase identi- 
fizierbar ist. 

Dieses stochastische Messmodell besteht aus einer analyti-- 
schen, deterministischen Messf unktion, welche die logarithm 
20 mierte, empfangene Signalf eldstarke in Abhangigkeit von Posi- 
tionskoordinaten beschreibt. Dartiber hinaus wird ein sto- 
chastischer TVnteil (Komponente) identif iziert, welcher Mo- 
de llunsiche r he it wie auch Messrauschen berticksichtigt . 

25 Ein solcher stochastischer Ansatz bei der Lokalisierxang be- 
rticksichtigt die Unsicherheiten bei der ermittelten Position 
durch Wahrscheinlichkeiten bzw. Wahrscheinlichkeitsdichte- 
f unktionen . 

30 Neben dem Messmodell fur die empfangbaren Kommunikationssig- 
nale wird die Gute der Lokalisierung dadurch weiter verbes- 
sert/ dass ein Benutzermodell mit eingebracht wird. Uber die- 
ses Benutzermodell kann Wissen tiber eine (raumliche und/oder 
zeitliche) Bewegung eines Benutzers mit berUcksichtigt wer- 

35 den. 
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Durch solche Benutzermodelle konnen maximale ^Schrittweiten^' 
von Benutzern in einem Zeitschritt, welcher tiber zwei zugeho- 
rige bzw. korrespondierende Signalmessungen (am Anfang \ind am 
Ende des Zeitschritts) definiert wird^ beschrSnkt ist. 

5 

Ein zusatzliche Verbesserung der Lokalisierungsgtite l^sst 
sich erreichen durch rekursive Positionsermittlung durch Kom- 
bination mehrfacher verschiedener nacheinanderfolgender Mes- 
sungen . 

10 

Der hier benutzte neue stochastische Ansatz bei der Lokali- 
sierung ftihrt zu nichtlinearen, mehrdimensionalen Messglei- 
chungen. Eine exakte LOsung dieses Hayes ianischen Filterprob- 
lems fiihrt zu komplizierten nicht-Gausschen Wahrscheinlich- 
15 keitsdichten^ welche die gesuchte Position mit ihren Unsi- 
cherheiten beschreiben. 

Diese Wahrscheinlichkeitsdichte (-funktionen) konnen durch 
Gaussche Mischdichten [19] r [20] angenahert werden. 

20 

Diese Annaherung wird hier durch einen neuentwickelten soge- 
nannten Prior Density Splitting Mixture Estimator'' (PDSME) 
durchge ftihrt/ der entwickelt wurde^ urn eine durch einen An- 
wender einstellbare GUte der Lokalisierung zu ermbglichen - 
25 im Gegensatz zu tiblichen, nichtlinearen^ bekannten Filter- 
techniken, wie einem Extended Kalman Filter (EKF) . 

Jedoch kann eine Angabe einer oberen Grenze fUr einen Rechen- 
aufwand durch BeschrSnkung der Anzahl von anwendbaren Gauss- 
30 Chen Mischkomponenten die Lokalisierungsgtite reduzieren. Dar- 
tlber hinaus wird durch die Representation von Dichten bei den 
Gausschen Mischungen eine rekursive Vorgehensweise sowohl bei 
den nichtlinearen ^Predictions Steps'' als auch bei den nicht- 
linearen "^Filtersteps" moglich. 

35 

Weiter wird im ersten Teil beschrieben werden, wie das sto- 
chastische Lokalisierungsproblem formuliert werden kann. Die 
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deterministischen Anteile und die stochastischen Anteile des 
Messmodells warden nachfolgend erklart werden. Die Anwendxing 
von Bayesianischen Filteralgorithmen, basierend auf eine 
gausschen Mischnaherung fur die exakten Wahrscheinlichkeits- 
5 dichten verwendet in dem neuen PDSME wird beschrieben werden. 

Im zweiten Teil werden Ergebnisse einer Lokalisierung eines 
DECT-Mobiltelefons in einem DECT-Netz dargestellt und die Gtl- 
te der Lokalisierung unter Verwendung des PDMSE mit der unter 
10 Verwendung eines bekannten EKF verglichen werden. 

Formulierung des Lokalisierungsproblems 

Die Lokalisierungsaufgabe fur ein mobile Kommunikationsein- 
15 richtung in einem Kommunikationsnetz kann in zwei Grundprob- 
leme eingeteilt werden: 

1) die Identif ikation des Messmodells mit einem determinis- 



neration phase; 

20 2) die Ermittlung der Wahrscheinlichkeitsdichte der gesuch- 
ten Position in einem Bayesianischen Filterstep in der 
Lokalisierungsphase . 

Das Messmodell 



beschreibt die Messung der logarithmierten empfangenen Sig- 
nalstarken 



des C~then Transmitters (Senders), C= 1/ . . . f N, als eine de- 
30 terministische nichtlineare Funktion hC(xk) der Positionsko- 
ordinaten xk=[xl,k, x2,k]fUr jeden der N verschiedenen Sen-- 
der, 

Der stochastische Anteil vk des Messmodells wird reprasen- 
35 tiert durch eine additive Unsicherheit . Additive Unsicherhei- 



tischen und einen stochastischen Anteil in der model ge- 



25 
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ten im Messmodell der logarithiaierten empfangenen Signalstar- 
ke entsprechen multiplikativen Unsicherheiten der Signalleis- 
tung welche auf einen Einfluss von iinbestimmten Dampfungen 
zurtlckzuftihren sind. 



Jeder Transmitter kann identif iziert werden tlber eine ihm ei- 
genen ID, welche durch die Basisstation wahrend einer Kommu- 
nikation tlbertragen wird. 

10 Dadurch wird das Lokalisationsproblem erheblich vereinfacht, 
weil jeder Messwert direkt der zugehSrigen nichtlinearen 
Messgleichiing zugewiesen werden kann. 

In der Lokalisierungsphase werden die nichtlinearen Filteral- 
15 gorithmen in jedem Zeitschritt k angewandt . Jedoch kOnnen oh- 
ne die additiven Unsicherheiten vk die impliziten, nichtline- 
aren Beziehungen zwischen dem Vektor lait den gemessenen loga- 
rithmierten empfangenen Signal star ken yk der Transmitter und 
den Positionen xk des Empfangers nur numerisch gelost werden, 
20 beispielweise durch eine Methode der kleinsten Fehlerquadra- 
te. BerUcksichtigt man die Unsicherheiten mit, so kSnnen die 
Messgleichungen verwendet werden, urn die gesuchte Position im 
Filterstep zu erneuern update"). 

25 Die weiter beschriebene Vorgehensweise bei der Lokalisierung 
basiert auf einer Naherung der exakten LOsung des bayesiani- 
schen Filtersteps durch Gaussche Mischdichten 



5 




mit L Komponenten, definiert durch nicht-negative Gewichte w, 
30 Erwartungswerte \x und Kovarianzen C. 



1st der Messwert y verfUgbar, so ergibt sich die exakte LO- 
sung des Bayesianischne Filtersteps aus: 
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mit der nichtlinearen Messfunktion h(x)/ den additiven Unsi- 
cherheitsdichtef unktiuonen f (v) \and eine Normalisieriingskon- 
stante c. 

5 A priori Wissen tibe die gesuchte Position wird reprSsentiert 
durch die Dichtefunktion f (x) . 

FUr jede neue Messung wird die vorangehende a posteriori 
Dichte f (x|y) interpretiert als neue a priori Dichte, so dass 
10 die gesuchte bzw. ermittelte Position rekursiv erneuert (up- 
date) werden kann. 

Analog zum Bayesianischen Filterstep wird ein Predikti- 
onsstep, welcher zur Beschreibung des nichtlinearen Benutzer- 
15 modells 

Xk+i = + (2) 

verwendet wird, ebenso zur Naherung der exakten Dichte durch 
Gaussche Mischungen durchgefilhrt • 

20 Die nichtlineare Funktion a(x) ist ein deterministisches Mo- 
dell einer Benutzerbewegung. Unsicherheiten werden wieder 
durch additives Rauschen w berUcksichtigt • 

Stochastische Modellierung 

25 

Das Messmodell wwird identifiziert tlber Messungen der loga- 
rithmierten Signalstarken fur jeden Transmitter an Gitter- 
punkten eines Gitters,, welchen die Lokalisierungsumgebung 
abdeckt. Die deterministische (h(x))und die stochastische 
30 Komponente (v) des Messmodells (1) werden vor der Lokalisie- 
rung identifiziert. 

Der deterministisch, analytische Anteil h(x) ist eine Mess- 
gleichung, welche die empfangene Signalleistung als eine 
35 Funktion der Positionskoordinaten beschreibt. Der stochasti- 
sche Anteil v ist ein Modell ftir die Unsicherheiten des de- 
terministischen Antteils. Diese Unsicherheiten setzen sich 
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zusammen aus raumlichen Unsicherheiten entsprechend Modell- 
fehlern und zeitliches Messrauschen. 

Deterministisches Messmodell der logarithiaierten empfangenen 
5 Signalstarken 

Fig.l zeigt, dass die Abnahme der logarithmierten, empfange- 
nen Signalstarke tlber eine Distanz von mehreren Metern anna- 
hernd linear ist. In anderen Dokumenten werden ahnliche An- 
10 nahmen beztiglich des Messmodells als ""linear-loss-iaodel" [24] 
Oder linear-scope-model'' [18] bezeichnet. 

Matliematisch lasst sich fur zweidimensionale Positionskoordi- 
naten x richtungsabhangige lineare Abnahme der logarithmier- 
15 ten empfangenen Signalstarken besclireiben durcli N unabhangige 
Me s s gl e i chungen 

= - + Ai; . (3) 

Die Parameter mC und PC der positiven semidef initen quadrati- 
schen Form und der additive Offset AC werden fUr jeden Trans- 
20 mitter C= 1/ .../ N identif iziert . 

Die 6N Parameter mCr PC und AC werden bestimmt durch eine Me- • 
tliode der kleinsten Felilerquadrate/ bei welcher die Abwei- 
cliung zwischen dem angenaherten Messmodell und der gemessenen 
25 logarithmierten Signalstarken auf den Gitterpunkte, gemessen 
waiirend der model generation phase, minimiert wird. Zur Re- 
duktion des zeitlichen Messrauschens wird der Mittelwert iiber 
mehrere Messungen an jedem Gitterpunkt zur Berechnung des 
Durchsctmitts der gemessenen Werte verwendet. 

30 

Stocliastische Modellierung der Unsictierlieiten 

Das stochastisclie Unsicherheitsmodell berticksichtigt sowolil 
die Abweicliung zwisclien den angenaiierten^ im obigen besclirie- 
35 benen Modell und der waliren Verteilung der logaritlimierten 
empfangenen Signalstarken,. gemessen auf den Gitterpunkten, 
und das zeitliclie Messrausclien . 
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Modellungewiliheiten 

Bei dieseia Lokalisierungsansatz wird die Abweichung zwischen 
5 der wahren logarithmischen Empf angssignalleistung und dem de- 
tearministischen Messungsmodell hlXk) durch eine Gaufische 
Rauschdichte mit Mittelwert ^i^^^ und Standardabweichung <t^^ 

fUr jeden Sender beschrieben. Diese GauBsche Ungewissheit ist 
eine Darstelliing des AnnSherungsf ehlers der deterministischen 

10 Komponente des Messungsmodells tlber die gesamte Lokalisie- 
riingsumgebung . Als erstes stellt es die mittlere Abweichung 
tlber die Lokalisierungsxjiagebung zwischen dem Modell und der 
wahren Empf angs lei stung aufgrund der unrichtigen Annahmen ftir 
den deterministischen Teil der Messungsgleichung dar. Die 

15 sich aus groBen Abweichungen zwischen der wahren Empf angs- 
signalleistung und dem Messungsmodell ergebenden Ungewisshei- 
ten sind in Fig. 1 dargestellt. Als zweites haben Messungen 
auch gezeigt, dass ortliche Abweichungen bestehen, die durch 
Reflektion, nichthomogene Ausbreitung der Funkwellen und In- 

20 terferenz verursacht sein konnten. Diese rSimlichen Variatio- 
nen der logarithmischen Empf angssignalleistung sind aus 
Fig. 2 ersichtlich, wo Messungen in einem Abstand von 2 cm 
gesammelt worden sind. 

25 Messrauschen 

Zusatzlich zu rSiamlichen Variationen der Empf angssignalleis- 
tung kann zeitliches Messrauschen durch Auswertung einer Rei- 
he von mehreren verschiedenen Messungen an einer festen Stel- 

30 le bestimmt werden. Dieses zeitliche Messrauschen ist eben- 

falls in Fig. 2 ftir drei verschiedene Messungen der logarith- 
mischen Empf angssignalleistung an jeder Messstelle darge- 
stellt. Dieses Rauschen wird wiederum durch eine GauBsche 
Dichte mit Mittelwert M'v^C Standardabweichung (T*^^ angena- 

35 hert. 
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Koiobination von Modellungewissheiten und Messrauschen 

Um ein einfaches Modell mit den beiden im obigen beschriebe- 
nen Ungewissheiten zu erhalten, wird weiterhin angenommen/ 
5 dass beide Ungewissheiten unabhSngig sind, Sie kOnnen daher 
durch eine einzige GauBsche Dichte /vpCC^'c) fUr jeden Sender 

? = 1,,..,N, modelliert werden, die definiert wird durch den 
Mittelwert 



Dieses Modell impliziert nicht nur unabhangige Ungewissheiten 
sondern es wird auch angenommen/ dass beide Ungewissheiten 
beschrieben werden kSnnen, ohne jegliche Positionsabhangig- 
keit in Betracht zu Ziehen. 

20 

Filteralgorithiaus 

In diesem Abschnitt wird eine kurze Ubersicht Uber den fUr 
die Lokalisierung von Funkkommunikationsvorrichtungen benutz- 
25 ten PDSME (Prior Density Splitting Mixture Estimator) gege- 
ben. Zusatzlich wird eine Anpassung der in Abschnitt 3 abge- 
leiteten Messgleichungen h{xjc) eingeftihrt/ um die Berechnung 
des PDSME-Schatzers zu vereinf achen, Weiterhin wird ein Pra- 
diktionsschritt fUr ein einfaches Benutzermodell dargestellt. 



10 



und die Standardabweichung 




15 



30 
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Die PDSME-Messungsaktualisier\mg 

Der in dieser Schrift dargestellte PDSME-Algorithmus fiir die 
Lokalisierung von Funkkommunikationsvorrichtungen beruht auf 
5 der Berechnung einer linearisierten Messungsaktualisierung 
ftlr Gaufische Mischungsdichten. Der Messiingsaktualisierungs- 
schritt dieses Filteralgorithmus ist in einem Blockschaltbild 
im oberen Teil der Fig. 3 dargestellt. Splitting (Aufspal- 
tung) beruht auf der Berechnung des Linearisierungsf ehlers 

0 



ftir jede Komponente der Gaufischen a-posteriori- 
Mischungsdichte, Dieses Kriterium ist dem Kullback-Leibler- 
15 Abstand [25] 

x.(r'.*')-/r(.«v).n(ffej^)^ 

zwischen der exakten a-posteriori~Dichte und ih- 

20 rer Annahervmg /a? iZklifkx) durch Ersetzen der nichtlinearen 
Messungsgleichung h(xjc) durch ihre Linearisierung am 
Mittelwert der iten Komponente der a-priori-Dichtfunktion 
sehr ahnlich. Durch Berechnen des Linearisierungsf ehlers (4) 
werden die am meisten zum Annaherungsf ehler der a-posteriori- 

25 Dichte beitragenden Gaufischen a-priori-Mischungskomponenten 

identif iziert . Zur Verringerung dieses Linearisierungsf ehlers 
werden diese Gaufischen Mischungskomponenten ersetzt^ indem 
sie in mehrere Mischungskomponenten mit kleineren Kovarianzen 
unter Verwendung von Auf spaltungsbibliotheken aufgespaltet 

30 werden, die off-line optimiert worden sind (siehe Fig. 4) . 



Nach dieser analytischen Neuabtastung" der a-priori~Dichte 
kann der Filterschritt durch eine Bank von EKF, Linearisie- 
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rung der Messungsgleichiing am Mittelwert jeder Komponente der 
Darstellung der Gaulischen Mischiing der a-priori-Dichte be- 
rechnet werden. 

5 Danach wird Redundanz durch einen Zusammenf tlhrungsschritt 
verringert/ indem mehrere Gaufische Mischiangskomponenten in 
eine einzige Gaufische mit vernachlassigbarem Annaherungsfeh- 
ler kombiniert werden. 

10 Anpassung des Messungsmodells ftir Vereinfachungen des Line- 
arisierungsf ehlerkriteriums 

Flir GauBsches Messungsrauschen kann der Linearisierungsf ehler 
(4) analytisch als lineare Kombination von Momenten der Dich- 

15 te Is (Sk\$kti) fUr polynomische Messungsgleichungen h^lxk) be- 
rechnet werden. In diesem Unterabschnitt wird daher eine An- 
passung der Messungsgleichung (3) abgeleitet. Nach einigen 
algebraischen Umwandlungen kann die quadrierte Messungsglei- 
chung. 

20 

CgfcC-A^)> \kk -ZStdpri + - £^)vc-vl) (6) 

als polynomische Funktion mit einem abgeSnderten "^Messwerf 
2jt^< und der transformierten Ungewissheit neu geschrieben 

25 werden. Aufgrund der nichtlinearen Transformation der Zu- 
fallsvariablen ist die Wahrscheinlichkeitsdichtefunktion 
/w,f nicht langer GauBisch. Im Lokalisierungsexperiment im 

Abschnitt 5 werden die exakten Momente erster und zweiter 
Ordnung von li^ berechnet, urn eine Gaufische AnnSherung von 
30 ftrXi^O bestiramen. Man beachte, dass die Momente von Sif 
vom Messwert S^m,^ abhSngig sind. Sie mtlssen daher ftir jede 
neue Messung der logarithmischen Empf angssignalleistung neu 
berechnet werden und sind nichit zeitlich unverSnderlich wie 
die im Abschnitt 3 beschriebenen Parameter jUu,^ xand au,^. 



35 
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Benutzermodellierung durch PDSME-Pradiktions-Schritt 

Ahnlich dem Filterschritt kann der PDSME auch auf die Berech- 
nung nichtlinearer Pradiktionsschritte angewandt werden. Ana- 
5 log zu nichtlinearen Filterschritten besteht die Berechnung 
eines angenaherten Pradiktionsschritts auch aus der Auswer- 
tung eines Linearisierungsf ehlers, einer Bank von lineari- 
sierten Pradiktionsschritten und der Verringerung der Anzahl 
von Gauflschen Mischungskomponenten in einem Zusammenftlhrungs- 
10 schritt (siehe unterer Teil der Fig, 3) . In dieser Arbeit 

wird nur ein lineares Benutze3cmodell in Betracht gezogen. Der 
Pradiktionsschritt kann daher analytisch berechnet werden, da 
die a-posteriori-Dichte durch eine GauBsche Mischungsdichte 
im Filterschritt angenahert worden ist. 

15 

Das PrSdiktionsmodell besteht aus der linearen Zustandsglei- 
chung 

20 

wobei der Mittelwert ^ des additiven Systemrauschens Wis- 
sen tiber mogliche Richtungen und mittJere 5chritticEngen der 
Bewegung des Benutzers darstellt. Die Kovarianzmatrix Co, von 
(0]c gibt eine Schatzung ftir die Verteilung der SchrittlSingen 
25 des Benutzers an. Ftir jede Komponente i = 1,...,L der GauJi- 
schen Mischungsdichte wird die vorhergesagte Gaufische Mi- 
s Chungs komponente dann beschrieben durch den Mittelwert 

30 

und die Kovarianz 



35 
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Die geschatzte Position kann dann als gewichtete Oberlagerung 



5 der Mittelwerte aller GauBschen Mischungskomponenten berech- 
net werden. 

Teil 2: 

10 Lokalisierungsexperiment 

In diesem Abschnitt wird ein reales Lokalisierxangsexperiment 
fUr die Validierung des beschriebenen Lokalisierimgsansatzes 
far DECT-Mobiltelefone dargestellt. In den Unterabschnitten 
15 wird die Oberlegene Leistung des PDSME im Vergleich ziiia EKF- 
Filter (Extended Kalman Filter) , einem weitverbreiteten stan- 
dardmafiigen Ansatz ftir die Zustandsschatzung von nichtlinea- 
ren Systemen, gezeigt. 

2 0 Auswertungss zenarium 

In diesem Lokalisierungsexperiment sind AT = 10 Sender in ei- 
nem Innenraum von annahernd 30 m x 30 m auf einem Stockwerk 
eines Gebaudes plaziert worden. Bei der Modellerzeugungsphase 
25 ist die logarithmische Empf angssignalleistiing jedes Senders 
auf einem Gitter mit 1 m Abstand zwischen Gitterpunkten ge- 
messen worden. Es sind W = 10 Messungsgleichungen h^iXk) , 
insgesamt bestehend aus 60 Parametern ftir die deterministi- 

schen Komponenten und 2 0 Parametern ftir die stochastischen 
30 Komponenten jj^.^ identif iziert worden. Aufler einigen wenigen 

Bereichen^ die durch hohe Dampfung der Funkwellen aufgr\ind 
von Stahlbetonwanden beeinflusst wurden^ ist das angenommene 
Modell eine zutreffende Annaherung der Verteilung der Emp- 
f angssignalleistung . 
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Die anfangliche Wahrscheinlichkeitsdichte der Position wird 
als GauBsche Dichte mit deia Anf angsmittelwert 



gewahlt^ d-h, es stehen beinahe keine Vorkenntnisse tiber die 
Position zu VerfUgung. 

Bei der Lokalisierungsphase wird eine Messung jedes empfang-- 
15 baren Senders entlang einer Linie wie in Fig. 4 dargestellt 

gesammelt/ mit einem Abstand von Ar = 0^5 m zwischen den wah- 
ren Messstellen Pi/ Pa/ • . • / -Pit • Die Positionskoordinaten der 
Messungen sind bestimmt worden, urn die Grundwirklichkeit mit 
den Schatzpositionen zu vergleichen. FUr die Lokalisierung 
20 des Empf^ngers sind weder im Lokalisierungsansatz mit dem 
PDSME noch mit dem EKF die genauen Messpositionen benutzt 
worden. Nachdem die Messungsaktualisierung fllr jeden empfang- 
baren Sender an einer festen Stelle Pi^i = 1/.../17 berechnet 
worden ist, ist ein PrSdiktionsschritt nach Unterabschnitt 
25 4.3 berechnet worden. Im vorliegenden Beispiel wird das ein- 
fache Benutzermodell definiert durch den Mittelwert 

Ho = £m 

30 und die Kovarianz 



und der Anf angskovarianz 



10 




d.h. es steht keine Kenntnis iXber eine bevorzugte Richtung 
35 der Bewegung des Benutzers zur Veffiigung. Durch die Kovari- 
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anzmatrix C^, wird nur die Verteilung der Schrittlangen des 
Benutzers dargestellt. 

EKF- Filter (Extended Kalman Filter) 

5 

Bei Verwendung des EKF- Filters wird die Messungsgleichung (5) 
am Mittelwert der a-priori-Dichtefunktion linearisiert . In 
der Fig. 5 sind die Ergebnisse des EKF fUr den ersten und den 
170sten Filterschritt dargestellt. Bei jedem Filterschritt 
10 wird die Messung der logarithmischen Empf angssignalleistung 
eines einzigen Senders zur Aktualisierung der geschatzten Po- 
sition benutzt. 

Es ist ersichtlich, dass es nach dem ersten Filterschritt ei- 
15 nen bedeutenden Schatzungsf ehler gibt. Die wahre Position des 
Empfangers liegt nicht innerhalb der StUtzwerte der geschatz- 
ten a-posteriori-Dichtefunktion, die durch ihre Profildar- 
stellung gezeigt wird, Weiterhin gibt es keine Kreuzungsstel- 
le zwischen der wahren Messungsgleichung und der geschatzten 
20 Dichte. Es ist daher nicht moglich, Datenvalidierungs- 

verfahren anzuwenden, urn f estzustellen, ob ein Messwert durch 
die geschatzte Dichte "^eriauterf werden kann. 

PDSME (Prior Density Splitting Mixture Estimator) 

25 

In der Fig. 6 werden die durch den PDSME-Algorithmus ge- 
schatztfen a-posteriori-Dichten auch fUr den ersten und 
170sten Filterschritt dargestellt. Die Annaherung der nicht- 
GauBschen a-posteriori-Dichte im ersten Filterschritt, die 
30 der numerisch berechneten optimalen Bayeschen Losung des Fil- 
terschritts sehr nahe liegt, ist eindeutig viel besser als im 
Fall des EKF- Filters. Es konnen daher nunmehr Datenvalidie- 
rungsverf ahren erfolgreich angewandt werden. 

35 In der Fig. 7 werden die Erwartungswerte der mit dem EKF und 
dem PDSME berechneten a-poster iori-Dichtefunktionen mit der 
wahren Position verglichen. Wenn man diese Figur mit Fig. 8 
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vergleicht, ist ersichtlich, dass fiir beinahe lineare Filter- 
probleme/ d.h. wenn die Kovarianz der Schatzposition im Ver- 
gleich mit der Nichtlinearitat der Messiingsgleichung nach 
mehreren Filterschritten klein ist, EKF und PDSME beinahe 
5 dieselben Ergebnisse bringen. In diesen Fallen benutzt der 

PDSME nur eine bescheidene Anzahl von GauBschen Mischungskom- 
ponenten wShrend fUr starke Nichtlinearitaten zu Beginn des 
Lokali si erungs experiments eine hOhere Anzahl von AnnSherungs- 
komponenten notwendig ist, vm den SchStzungsf ehler zu verrin- 
10 gern. Die Uberlegene Leistung des PDSME-Filters im Vergleich 
mit dem EKF wird auch durch den Durchschnitts-Schatzf ehler 



15 iiber die Nf == 170 Filterschritte gezeigt, der 3,30 m ftir den 
EKF und 1,22 m ftlr den PDSME betrSgt. 

Schlussfolgerungen 

20 In dieser Arbeit ist ein stochastischer Ansatz fUr die Loka- 
lisierung von Funkkoramunikationsvorrichtungen geboten worden, 
der auf dem Messen der logarithmischen Signalleistung der 
empfangbaren Sender durch eine Mobilstation beruht. Ftir jeden 
Sender ist ein stochastisches Messungsmodell identif iziert 

25 worden, das aus einer deterministischen und einer stochasti- 
schen Komponente besteht. Dieses Messungsmodell ist far die 
Schatzung der Position des EmpfSngers durch einen neuartigen 
GauBschen Mischungsschatzer benutzt worden, der auf dem Auf- 
spalten der a-priori-Dichte entsprechend einem Linearisie- 

30 rungsfehlerkriterium beruht- Dieses Kriterium ist dem Kull- 

back-Leibler-Abstand zwischen der wahren und der angenaherten 
a-posteriori-Dichte sehr ahnlich, die durch eine Linearisie- 
rung der Messungsgleichung berechnet wird, Durch Anwendung 
dieses Schatzungsverf ahrens auf die Lokalisierung von DECT- 

35 Mobiltelefonen konnen bedeutende Verbesserungen der Schat- 
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zungsgiite erreicht werden, wenn statt der standardmaBigen An- 
satze wie dem EKF der PDSME benutzt wird. Weitere Verbesse- 
rung der Lokalisierungsgiite kann durch Identif izieren besser 
deteinninistischer Messungsmodelle iind genauere Keimzeichnun- 
gen des Messungsrauschens erreicht werden, die keine Unabhan- 
gigkeit zwischen den in dieser Arbeit erwShnten verschiedenen 
Ungewissheiten annehmen. 
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Patentanspruche 

l.Verfahren zur Ermittlung einer ausgewahlten Position einer 
mobilen Koiniaunikationseinrichtung in einem Kommunikationsnetz 
mit mindestens einer ersten Basisstation, eingerichtet zu ei- 
ner ersten Koininunikation mit der mobilen Kommunikationsein- 
richtung mittels erster Kommunikationssignale/ 
unter Verwendung eines nichtlinearen Kommunikationsmodells 
mit einem deterministischen Anteil und einem stochastischen 
Anteil/ 

- welcher deterministische Anteil eine Abhangigkeit zwi- 
schen Koiraaunikationssignale einer Basisstation und Posi- 
tionen der mobilen Kommunikationseinrichtung beschreibt 
und 

- welcher stochastische Anteil eine Unsicherheit des de- 
terministischen Anteils beschreibt, 

- bei dem unter Verwendung des zu der ausgewahlten Position 
zugehSrigen ersten Kommunikationssignals der ersten Basis- 
station und dem nicht-linearen Kommunikationsmodell die 
ausgewahlte Position ermittelt wird/ 

- wobei bei der Ermittlung eine nichtlineare Bayesianische 
Filtertechnik angewandt wird, wodurch ein mSgliches Auf- 
enthaltsgebiet der ausgewahlten Position bestimmt wird^ 
aus welchem die ausgewahlte Position ermittelt wird. 

2.Verfahren nach Anspruch 1, 

bei dem' mindestens eine zweite Basisstation, eingerichtet zu 
einer zweiten Kommunikation mit der mobilen Kommunikations- 
einrichtung mittels zweiter Kommunikationssignale, verwendet 
wird, wobei unter Verwendung des zu der ausgewShlten Position 
zugehorigen ersten Kommunikationssignals und/oder des zu der 
ausgewahlten Position zugehSrigen zweiten Kommunikationssig- 
nals und dem nicht-linearen Kommunikationsmodells die ausge- 
wahlte Position ermittelt wird. 
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3.Verfahren nach Anspruch 1 oder 2, 

bei dem unter Verwendung des ersten und/oder des zweiten Kom- 
munikationssignals ein und/oder mehrere entf ernxmgsabhangige 
Parameter ermittelt warden bzw. wird, welche bzw. welcher von 
5 einer Entfernung der mobilen Kommxanikationseinrichtiing zu ei- 
ner von den Basisstationen abhangen bzw, abhangt und unter 
Verwendung dessen oder derer das mOgliche Auf enthaltsgebiet 
ermittelt werden bzw. wird. 

10 4.Verfahren nach einem der vorangehenden AnsprUcher 

bei dem die Unsicherheit des deterministischen T^nteils eine 
Unsicherheit der Kommunikationssignale/ insbesondere ein 
Messrauschen, und/oder eine Unsicherheit der libhangigkeit/ 
insbesondere eine Unsicherheit des Kommunikationsmodells, 

15 sind bzw. ist. 

S.Verfahren nach einem der vorangehenden Ansprtlche/ 
bei dem das zu der ausgewahlten Position zugehorige erste 
Kommunikationssignal zeitabhangig ist und/oder fUr einen 
20 Zeitpunkt k gemessen wird. 

S.Verfahren nach einem der vorangehenden Ansprtlche^ 

bei dem die nichtlineare Bayesianische Filtertechnik iterativ 

angewandt wird. 

25 

T.Verfahren nach einem der vorangehenden Ansprtlche, 
bei dem als die nichtlineare Bayesianische Filtertechnik ein 
Gausscher Mischf ilter-Algorithmus oder ein Kalman Filter, 
insbesondere ein Extended Kalman Filter, verwendet wird. 

30 

S.Verfahren nach Anspruch 7, 

bei dem der Gaussche Mischf ilter-Algorithmus ein Prior Densi- 
ty Splitting Mixture Estimator (PDSME) ist. 

35 9.Verfahren nach einem der vorangehenden Ansprtiche, 

bei dem ein Benutzermodell verwendet wird, welches eine Bewe- 
gung der mobilen Kommunikationseinrichtung beschreibt. 
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10. Verfahren nach Anspruch 9, 

bei dem unter Verwendung des Benutzermodells die Bewegung der 
mobilen Kommunikationseinrichtung fUr einen Zeitschritt be- 
grenzt wird, 

5 

11. Verfahren nach Anspruch 9 oder 10, 

bei dem unter Verwendung des Benutzermodells eine Positions- 
schatzung fUr die ausgewahlte Position ermittelt wird. 

10 12. Verfahren nach Anspruch 11, 

bei dem die Positionsschatzung unter Anwendung der nichtline- 
aren Bayesianischen Filtertechnik ermittelt wird. 

13. Verfahren nach Anspruch 11 oder 12, 
15 bei dem die Positionsschatzung als ein Startwert fiir die Er- 
mittlung des mogliches Auf enthaltsgebiet der ausgewahlten Po- 
sition durch die Anwendung der nichtlinearen Bayesianischen 
Filtertechnik verwendet wird. 

20 14. Verfahren nach einem der vorangehenden Ansprilche, 

bei dem das Kommunikationsnetz ein WLAN-, GSM- oder DECT-Netz 
ist- 

15. Verfahren nach einem der vorangehenden Ansprtlche, 
25 bei dem als ausgewahlte Position der mobilen Kommunikations- 
einrichtung ein Kennwert des maglichen Auf enthaltsgebiets, 
insbesondere ein Schwerpunkt oder ein Erwartungswert des m5g- 
lichen Aufenthaltsgebiets, verwendet wird. 

30 16. Verfahren nach einem der vorangehenden Ansprtiche/ 

eingesetzt zur Lokalisierung eines Mobiltelefons in einem di- 
gitalen, zellularen Mobil funknetz, insbesondere in einem GSM- 
Netz, 

wobei die mobile Kommunikationseinrichtung das Mobiltelef on, 
35 die erste Basisstation eine gesprachsfiihrende Basisstation in 
dem Mobilfunknetz ist. 
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17.Anordnung zur Ermittlung einer ausgewahlten Position einer 
mobilen Kommunikationseinrichtung in einem Kommunikationsnetz 
mit mindestens einer ersten Basisstation, eingerichtet zu ei- 
ner ersten Kommunikation mit der mobilen Kommunikationsein- 
5 richtung mittels erster Kommunikationssignale, 

unter Verwendxing eines nichtlinearen Kommunikationsmodells 
mit einem deterministischen Anteil und einem stochastischen 
Anteil^ 

- welcher deterministische Anteil eine Abhangigkeit zwi- 
10 schen Kommunikationssignale einer Basisstation und Posi- 

tionen der mobilen Koiranunikationseinrichtung beschreibt 
und 

- welcher stochastische Anteil eine Unsicherheit des de- 
terministischen Anteils beschreibt, 

15 - mit einer Positionsermittlungseinheit/ die derart einge- 
richtet ist, dass unter Verwendxing des zu der ausgewahlten 
Position zugehorigen ersten Kommunikat ions signals der ers- 
ten Basisstation und dem nicht-linearen Kommunikationsmo- 
dell die ausgewShlte Position ermittelbar ist, 

20 - wobei bei der Ermittlung eine nichtlineare Bayesianische 
Filtertechnik angewandt wird, wodurch ein mogliches Auf- 
enthaltsgebiet der ausgewahlten Position bestimmt wird, 
aus welchem die ausgewahlte Position ermittelt wird. 

25 IS.Computerprogramm mit Programmcode-Mitteln/ um alle Schrit- 
te gemafi Anspruch 1 durchzuftihren, wenn das Programm auf ei- 
nem Computer ausgeftihrt wird. 

19. Computerprogramm mit Programmcode-Mitteln gemafi Anspruch 
30 18, die auf einem computerlesbaren Datentrager gespeichert 

sind. 

20. Computerprogramm-Produkt mit auf einem maschinenlesbaren 
Trager gespeicherten Programmcode-Mitteln, um alle Schritte 

35 gemafi Anspruch 1 durchzufuhren, wenn das Programm auf einem 
Computer ausgeftihrt wird. 
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